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The chiral extension of Quantum Chromodynamics (XQCD) adds to the standard lattice action explicit pseu- 
doscalar meson fields for the chiral condensates. With this action, it is feasible to do simulations at the chiral 
limit with zero mass Goldstone modes. We review the arguments for why this is expected to be in the same 
universality class as the traditional action. We present preliminary results on convergence of XQCD for naive 
fermions and on the methodology for introducing counter terms to restore chiral symmetry for Wilson fermions. 
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1. INTRODUCTION 

In spite of the impressive progress in extract- 
ing useful physics from the lattice formulation of 
Quantum Chromodynamics, in almost all cases 
the results rely on the quenched approximation. 
The additional cost of simulating full QCD with 
internal fermion loops forces one to unrealisti- 
cally small lattices and heavy masses for the light 
quarks. Even with the advent of Teraflops scale 
computing, new methods will be needed to avoid 
these limitations. Improved methods should be 
sought (i) to allow accurate results on smaller 
lattices (i.e. larger lattice spacings at a fixed vol- 
ume) and (ii) to extrapolate reliably to the chiral 
regime of light quarks. 

Here we discuss a new lattice action (XQCD), 
which is intended to help with the second problem 
(ii) by extending lattice QCD to include explicit 
fields for the low mass Goldstone modes: the pion 
and eta. To appreciate the concept behind this 
extension, it is useful to consider it as a part of a 
more general approach that includes recent efforts 
to address the first issue (i). 

As originally suggested by Wilson in his renor- 
malization group (RG) approach, higher dimen- 
sional operators should be added to a lattice ac- 
tion to suppress scaling violations due to the lat- 
tice spacing, a, or the finite momentum cut-off. 
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A = Tr/a. Recently this program has gained new 
momentum because of the realization that even 
a few higher dimensional operators can substan- 
tially reduce the lattice artifacts. For example, 
the clover action |^ adds to the five dimensional 
Wilson term, A~^'!/;D^i?^'0, the only other five 
dimensional operator. 
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With the proper redefinition of the fermion 
fields, this allows one to remove the 0{a) and 
O(a5'o"^og"(a)) lattice artifacts in the quark 
propagator. 

At dimension six, the pure gauge operators. 



A- 



'Tr[D^F^^,D^F^^,], 
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when expressed in terms of the appropriate Wil- 
son loops, lead to improved gauge actions. From 
the study of toy models, there is also optimism 
that simple approximations to a classically scale 
invariant (or "perfect" ) action can even more dra- 
matically reduce lattice artifacts, for asymptoti- 
cally free theories 

At first sight our chiral extension of QCD ap- 
pears to be a departure from this trend in renor- 
malization group improved actions. However, if 
we simultaneously block fermionic and gauge de- 
grees of freedom, we must bring in higher polyno- 
mials in the fermion fields as well. For example, 
in addition to several new gluon quark-bilinears 
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similar to Eq. (|l|), the remaining six dimension 
terms are local four-fermion operators. With two 
flavors, enforcing U{2) x U{2) chiral-flavor invari- 
ance, there are only 4 possible four-fermion oper- 
ators: the pion-eta operator, Oe{x), 



A" 



[(V;V')' + (V^i75V')'+ {^pT^)^ + {^ni5T^j)% (3) 



and three others, A~^{^pjfj,ipy , A~^{'ipi^5jf^^)'^ , 
and A~^[ {ip^^Tij})'^ + (V'*757/j'?V')^ which are of 
less interest to us at present. 

In its present form, XQCD amounts to a gen- 
eralized lattice action resulting from bosonizing 
the first of the four-fermion operators, (||). There 
are may other ways to motivate our choice. For 
example in the standard discussions of the ori- 
gins of chiral symmetry breaking, one arrives at 
this same four-fermi operator (||) by the most at- 
tractive channel argument for multiple gluon ex- 
change graphs. By introducing instanton effects, 
one is led to the simpler four-fermion term of the 
Nambu-Jona-Lasino (NJL) model. 
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which is a linear combination of Oq{x) and the 
'tHooft determinant. The explicitly breaking chi- 
ral U(l) symmetry in NJL model removes the 
need for an eta Goldstone mode. Although we 
are studying XQCD for both cases, Eqs (y) and 
(0), it is probably better to avoid introducing ex- 
plicit chiral U(l) breaking at this point, since on 
the lattice the anomaly is properly accounted for 
by the doublers as their masses go to infinity with 
the cut-off. 

Eventually, we may consider the complete set 
of six dimensional operators using renormaliza- 
tion group methods to fix their coefhcients. In- 
deed, several years ago ^, we began with the 
more ambitious goal of deducing XQCD by RG 
transformations of QCD on a fine lattice. How- 
ever, it is also useful to investigate the immediate 
consequence of the simplest possible chiral exten- 
sions of QCD in the same spirit as the limited in- 
vestigations of separate improved actions for the 
pure gluon and the Wilson fermion sectors. 



2. XQCD EFFECTIVE ACTION 

The standard lattice Lagrangian as formulated 
by Wilson is 

I^QCD - ^i?eTr[l - C/^,,(x)] -f 7/iM(C/)^ (5) 
with fermionic matrix. 
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^^[2^.t. 



rriQ-il^^ii^ (6) 



To simplify the present discussion of chiral invari- 
ance, we drop the Wilson term (r — 0) and ignore 
the doubling problem. 

The Wilson QCD action has a UV fixed point 
at go = 0- By tuning the bare gauge coupling to 
the UV fixed point 0, the lattice correlation 
length, ^, will increase and eventually we enter 
the scaling region where the lattice theory de- 
scribes a continuum QCD theory with renormal- 
ized coupling g. One may add our four-fermion 
operator (H) to the bare QCD action with an ar- 
bitrary dimensionless coupling Go, 

GMx) ■ (7) 

According to the RG theory, up to 0{^^^ In'' ^) 
terms, we can always match the renormalized the- 
ory of the actions in Eq. (||) and Eq. (|^) such that 
the effect of the four-fermion operator can be ab- 
sorbed in a suitable choice of the bare gauge cou- 
pling. In this sense the higher dimensional oper- 
ators are "irrelevant" : the two actions belong to 
the same universality class and their continuum 
limits are equivalent. 

By a standard procedure, the four-fermion op- 
erator can be "bosonized" 
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introducing a complex Lagrange multiplier field 
^x — 4>o + if ■ (f), where the tilde denotes the 75 
projections. 



1 + 75 



1 
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(9) 
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and j/o = V2G0. It is very important to note that 
the "mass" parameter for the scalar field is 0(A), 
diverging with the cut-off. 
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Now we may imagine doing additional RG 
blocking transformations on the new action. All 
operators that are absent in Eq. but allowed 
under the chiral and gauge symmetries will be 
generated. In particular, both a kinetic term and 
a chirally invariant quartic term will be generated 
for the <i> field, 

^igczj = ^QCD + Ci + yo i^Ai'. , (10) 
where 

+ y(rr[$t$,]). (11) 

It is this generalized scalar extension that we pro- 
pose as a candidate action for XQCD. The scalar 
field propagates with a mass at the cut-off so 
only short range polynomial quark interactions 
are generated by its exchange. 

On the lattice, it is also convenient to re-write 
the scalar field $ as V^^^S, where S is an U(2) 
matrix field. It is well known in lattice Higgs 
theories that one can tunc the potential so that 
the radial mode is frozen, 

$^S = e*^/^'' e^^-^/^-, (12) 

without changing the universality class. We will 
do this to reduce the number of new degrees of 
freedom. We are studying two cases: U(2) ex- 
tended QCD with the eta field and SU(2) ex- 
tended QCD without the eta field. 

2.1. Induced Chiral Symmetry Breaking. 

XQCD is constructed from two chiral models: 
lattice QCD and the S model coupled through 
the Yukawa term. Consequently for ^ there 
is a single U(2) chiral symmetry and there must 
be two phases - a symmetric phase and a bro- 
ken phase possessing the Goldstone bosons with 
the quantum numbers of the tt and 77. In the 
broken phase the physical Goldstone modes are 
mixtures of the "elementary" E and "composite" 
("075 V': ■0*75 ^"0) chiral fields. 

Here we give a qualitative survey of the phase 
diagram in the (/?, k, j/q) parameter space. Since 
EJ,S]j; = 1,2/0 has units of mass. For definiteness, 
we restrict the discussion to the SU(2) version 



of XQCD (dei(E) = 1) without the r] Goldstone 
mode. 

For j/o = 0, the scalar field sector, which is 
decoupled from the QCD, is the 0(4) symmetric 
model with a critical point at Kcr- The broken 
phase {k > Kcr) has a non- vanishing order pa- 
rameter (T = iTr[E]. In the symmetric phase 
{ K < Kcr), the mass gap (or sigma mass squared) 
is given by the difference, 

a^m^ = 1/k— I /Kcr- (13) 

We can anticipate that this is the appropriate re- 
gion for regaining the correct continuum limit of 
QCD, since we expect on the basis of our ear- 
lier discussion that the mass squared parameter 
should be positive (not tachyonic) and of the or- 
der of the cut-off ( rric ^ A — n/a ). 

At 2/0 = we know that chiral symmetry in 
the QCD sector is broken (< 'tptp >^ 0) for all 
values of /3. Therefore, as indicated in Fig. |l| by 
the horizontal dotted line, we have two phases in 
the 2/0 = plane: 

< a > = 0, < -00 >^ 0, K < Kcr , 

< a > ^ 0, < ■0?/' >7^ 0, K> Kcr ■ 

In the 50 = plane, Eq. (8) becomes the 
SU{2) X SU{2) symmetric Higgs- Yukawa model, 
which has been studied extensively For our 
purpose, we are only concerned with the critical 
surface in the weak 2/0 region (see Fig. |^). The 
symmetric and broken phase is separated by a 
second order phase transition line, Kcrivo), con- 
nected to K = Kcr Btt yo — 0. 

When we couple the two models (2/0 7^ 0), it is 
important to determine if the disorder of the E 
sector can defeat the chiral breaking mechanism 
of the standard QCD action. There are no numer- 
ical results available in the region where (/3, k, 2/0) 
are all finite. However an expansion in powers 
of 2/0 can be easily performed, since this is the 
only coupling between the QCD and the E sec- 
tors: SxQCD = SqcD + S^+yo J2x 'ipx^xlpx- 

To first nontrivial order in 2/0 : the order pa- 
rameters < a > and < 'ip^p > for chiral symmetry 
breaking are 

(a) = ((7)e ~ yo Xa {■4>^)qcd + 0{yl) , (14) 
(■00') = (0'0')qcd - yo(o-)s Xviv + '^(yo)- (15) 
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Figure 1. Phase diagram for XQCD. 



The VEV's are ensemble averages in the QCD 
and S sectors respectively, 



(O) 



QCD 



= [I/Zqcd) [dUdi>di>]0 



-Sqcd 



(1/Zs) /[dqoe 



(Oh 



and and Xc are the corresponding suscepti- 
bilities. 

Since we know that the QCD chiral symmetry 
is always broken {{'ipip)QCD 7^ ) for < /3 < cxd , 
it is clear from Eqs. (^) and (|l|) that XQCD 
must be in the broken phase < cr >7^ 0, < ip^p >^ 
at any /3 and k values in the small yo region. 
Intuitively, this result is easy to understand. A 
paramagnetic material will have nonzero mag- 
netization in an external magnetic field. For 
XQCD, the QCD sector acts as an external sym- 
metry breaking source to the scalar sector. Con- 
sequently, we are always guaranteed to be in the 
correct chiral phase for continuum QCD. To reach 
the continuum limit, we must approach the criti- 
cal surface at = without simultaneously tun- 
ing to the second order line at Hcrivo)- This addi- 



tional "fine tuning" would introduce an unwanted 
finite mass scale for the Higgs field. 

We have also explored the boundaries of the 
phase volume at the strong gauge coupling plane 
f3 = and at the strongly disordered chiral plane 
K = to check that there is no sign of chiral sym- 
metry restoration. We have carried out a per- 
turbative study for the weak coupling limit of 
XQCD As a result, we are confident 

that the only critical surfaces lie in the weak cou- 
pling plane at ^ g"^ — Q 'as depicted in Fig. |l|. 

2.2. Toy Model Results 

XQCD resembles earlier chiral quark mod- 
els, such as the Georgi-Manohar (CM) model, 
in which both quarks and explicit pions have 
been introduced. However, unlike these models, 
XQCD is not meant to be a weak coupling effec- 
tive theory. We do not assume new mass scales 
as implied by the condition, K^sb « ^qcd, in- 
troduced in the CM model. In XQCD, all "spuri- 
ous" degrees of freedom must have masses at the 
cut-off scale. To understand how these properties 
might come about in a non-perturbative context, 
in a recent paper , we have investigated several 
toy models in the large N limit. 

A natural toy model for XQCD is the Extended 
Nambu-Jona-Lasino model (XNJL), which cou- 
ples the sigma model to a NJL model via a 
chirally invariant Yukawa term, yoip^a + i^tsT^ • 
Tr)ip. Similar to our expectation for XQCD, the 
physical sigma resonance and pion in this ex- 
tended model occur as a linear combination of 
elementary and composite operators: (cr, ■?) and 
(jpip , 'ipij^Tip) , respectively. 

We have demonstrated that the pseudoscalar 
mass spectrum is. 



mi = 



m 



2 ^ f.2y0i<tip_> <o_^, 



(Here, is the four-fermion coupling of the NJL 
model.) There is only one massless pion and the 
extra pion, tt', resides at the cut-off scale. Similar 



results hold for the scalars. One mass, 



which 



is proportional to the VEV, remains in the phys- 
ical spectrum, while the other mass, m^' ~ fn-^', 
goes to infinity with the cut-off. Finally, we have 
verified that as the Yukawa coupling, yo , is tuned 
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to zero the spurious states remain at the cut- 
off. Just as we conjecture for XQCD, there is 
a smooth decouphng hmit, where only the physi- 
cal pion remains in the spectrum. In the limit of 
infinite cut-off XNJL and NJL models (at least 
at large N) are identical. 

2.3. Constituent Mass Extrapolation. 

We have started to do numerical simulations 
of XQCD with naive fermions using the 4> algo- 
rithm weighting the Trace-Log term appro- 
priately for Hf — 2. The simulations were per- 
formed on a 32-node CM5 using the CDPEAC 
routines of the BU-MIT collaboration code. For 
this preliminary investigation, the lattice size was 
8^ at /9 = 5.6 and k = 0.3. For QCD, we varied 
the fermion mass toq from 0.005 to 0.1, while for 
XQCD, we kept mo = 0.0 and varied yo from 0.01 
to 1.0. All runs started from a cold Higgs field and 
a quenched gauge configuration at /3 = 5.8. Af- 
ter performing 100 trajectories each with 30 MD 
steps, we measured the number of iterations for 
the conjugate gradient (CG) algorithm to con- 
verge to a residue of 10~^^. The results are plot- 
ted in Fig. |. 

One concludes that we can simulate naive 
fermions exactly at the chiral point, mg = 0, 
without any convergence problem. Furthermore 
the convergence rate of CG for XQCD depends 
on the effective constituent mass yo < cr > oi 
the fermions and it exhibits the same convergence 
rate as QCD compared atrriQ ~ yo < a >. Larger 
lattices with a variety of values for k,?/o) are 
under investigation to better estimate the relative 
efficiency between XQCD and QCD simulations. 

The difficulty encountered with the standard 
QCD actions is the necessity to extrapolate from 
large bare quarks masses for u and d quarks 
(~ 50-200 MeV) to the physical masses ( 5-10 
MeV). In XQCD, we may consider an alterna- 
tive extrapolation in the constituent quark mass 
by taking the limit yo ^ 0. Obviously in this 
limit, XQCD becomes the standard lattice the- 
ory. In fact this approach is not so silly, if the 
limit is smooth and none of the spurious state 
are re-introduced into the low mass spectrum as 

On the other hand, it should not be necessary 
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Figure 2. Comparison of the convergence rate for 
the conjugate gradient algorithm in QCD versus 
XQCD. 



to extrapolate to yo — 0. If we are in the same 
universality class, the fact that we regain con- 
tinuum QCD as the lattice spacing goes to zero 
means that yo is an irrelevant parameter and the 
physical spectrum should be only weakly depen- 
dent on the value of yo over a range of small val- 
ues of yo- We are investigating numerically this 
range, but on intuitive grounds any value of yo 
which contributes a relatively small part of the 
300MeV constituent mass should be safe. 

3. DOUBLER PROBLEM 

At first, it appears that a Wilson- Yukawa sec- 
ond order derivative term would be the ideal 
choice for removing doublers, because it preserves 
chiral invariance. At go = 0, this leads to the 
Smit-Swift model which was studied extensively 
in the last few years. There it was found that 
a Wilson- Yukawa term either can not lift the 
doublers or removes all fermions from the spec- 
trum llTtl. 
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However, there is one difference wtiich ap- 
peared at first to allow a new solution. In our 
application, the VEV associated with the Higgs 
field can be taken to the cut-off. We have stud- 
ied this scenario in the context of the XNJL toy 
model and found again that the spectrum looks 
fine, with one light pseudoscalar meson and both 
the doublers and the spurious scalar states at the 
cut-off. But there continues to be a fatal disease 
that looks to be very general. Namely in spite 
of the separation in mass scales, the value of F^^ 
for the Goldstone mode remains at the cut-off. 
Consequently we are not sanguine about this ap- 
proach. 

We are left with the options to use either the 
original Wilson term or staggered fermions. 
The obvious problem with both choices is their 
explicit breaking of chiral U(2) symmetry. Both 
options require counter terms to restore (in the 
continuum limit) full chiral symmetry. 

Before considering Wilson fermions in more de- 
tail, we remark that staggered fermions introduce 
very similar issues. One can couple the S fields 
symmetrically to pairs of staggered flavors so that 
by the usual trick of taking the square root of the 
fermionic determinant, one is left with an SU(2) 
flavor form of XQCD. One nice feature of stag- 
gered fermions is that the U(l) center of the chi- 
ral SU(2) group is preserved so that vector Ward- 
Takahashi identities (WTI's) alone can be used to 
restore the entire SU{2) x SU{2) current algebra. 
Nonetheless, we feel that the price one pays for 
the flavor symmetry breaking is still high and we 
prefer to use Wilson fermions. 

4. COUNTER TERMS 

Since the Wilson term breaks chiral invariance, 
the Ward-Takahashi identities (WTI's) associ- 
ated with chiral symmetry are no longer satisfied. 
One of these chiral Ward identities is the vanish- 
ing of the pion two point function at zero external 
momentum G~^(0) = 0. For the standard Wilson 
formulation, it can be shown that to restore chi- 
ral symmetry in the continuum limit this is the 
only condition that must be met. Thus the non- 
perturbative mass counter term is introduced by 
tuning the bare quark mass mo = Co + niq, 



where mo = Cq for a zero mass pion. (This is 
usually expressed in terms of the Wilson hopping 

parameter, K = Kcrit + 0{a).) 

4.1. Axial Vector Current 

To restore chiral symmetry to XQCD, we must 
look at additional WTI's. The full set of WTI's is 
found by local chiral rotation for the integration 
variables in the path integral. For example, SU(2) 
chiral invariance leads to 

< O Af.A^^ix) >=< O V^75T^(mo + xa)^ > (16) 
where 

+ \nTr[{Ta,^l}^^+^]-h.c. (17) 

For simplicity we have taken O to be chiral invari- 
ant. The xa term on the right comes from the 
variation of the Wilson term. To impose current 
conservation, we must introduce counter terms to 
cancel the singular and constant contributions of 
Xa in the limit of zero lattice spacing. 

As one may easily verify in lowest order per- 
turbation theory, the set of local counter terms 
required to restore chiral symmetry is: 

4 

VcT = ^C„rr[(E)"]-f C'5rr[(S,+^-E,)2] 

+ h.c. + Alpil; + B tpi:''tp. (18) 

There are additional counter terms for the wave 
functions and the axial current, but these play 
a somewhat different role since they do not have 
to be introduced into the bare action prior to the 
simulations. They are only needed at the analysis 
stage. 

4.2. Non-Perturbative Counter Terms 

To fix the symmetry restoring counter terms, it 
is sufficient to consider two matrix elements ( the 
vacuum and the single quark matrix elements) 
of the WTI dH) with the Ex field treated as an 
external source. Since the singular contributions 
only occur for orders up to j/qi it is convenient 
to expand for small yo. By expanding in j/Qj the 
path integral factorizes between the QCD and the 
E sectors (see Section 2.1 ). 
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For example for the vacuum matrix element, 
one just computes term by term (up to Uq) the 
vacuum expectation values of 

Trh^r^XA^^ + + Et) + lyoj,{^ - St)] 

in QCD at a fixed renormalized quark mass, 
iriq = /i, computing the necessary traces by us- 
ing the standard method of averaging over Gaus- 
sian pseudo fermionic sources. Each coefficient, 
Cniaya,gQ,a^,a), is then expressed as a low or- 
der polynomial in yo, so that once it has been de- 
termined, XQCD is defined for a range of values 
of yo without re-computing the counter terms. 



Table 1 

One loop fermionic integrals. 



/„ 


m 


m = 1 


m = 2 


m = 3 


m = 4 


n= 


= 1 


0.2347 


0.8486 


3.5421 


16.318 


n= 


= 2 


0.0259 


0.0597 


0.1925 


0.7366 


n= 


= 3 


(log-div) 


0.0086 


0.0165 


0.0470 


n= 


= 4 


(lin-div) 


(log-div) 


0.0031 


0.0050 



We are in the process of computing the counter 
terms non-perturbatively and comparing them 
with zeroth and first order expansions in cjq. To 
zeroth order, the single fermion loop contributes 
to C„, n = 1, .., 5: 



C: 



1 



a" a 
2 = -^yll2.2 + 2y^{h.2-hA) + 0{9l), 



C3 = ^3,3 + O(go'), 

Cs = vllz + Oigl) 

The single loop integrals are, 

'■'^ d'^p [rw{p) + anY 



, (27r)4 [g{p)]n 



(19) 



with g{p) = + {'>'w{p) + afi)"^ and 

w{p) — J2fj.i^ ~ cos{p^j)). The values at /i = 



and r = 1 are given in Table |. Since the Wil- 
son term vanishes at low momentum, no IR di- 
vergences occur in the chiral breaking terms to 
this order. The expression for the wave function 
renormalization integral, Iz, is a little more com- 
plex but it has the finite value Iz = 0.05473. The 
first contributions to A and B are 0{g1). 

For U(2) extended QCD, we also need to im- 
pose the WTI for the U(l) axial current. Due to 
the anomaly, this case is a little more intricate so 
we postpone it to a future publication. 

5. CONCLUSIONS 

The chiral extension of QCD, provides an im- 
proved action in which one is able to simulate 
near to or even at the chiral limit. Thus one 
can study the chiral vacuum at the physical mass 
of the pion. To avoid doublers, in the Wilson 
scheme, counter terms must be introduced, but 
they can be computed numerically in standard 
QCD at the decouphng point yo = 0. By studying 
the dependence of XQCD on the Yukawa coupling 
one can test its insensitivity to y^, an important 
signal of universality. Also extrapolations in the 
constituent mass of the quark, afford an al- 
ternative to the current practice of extrapolating 
in the bare quark mass, toq. We believe that be- 
ing able to invert the XQCD quark propagator 
in the chiral limit is an important advantage of 
the chiral extension. Simulations are underway 
to determine if this is true. 

REFERENCES 

1. G. Heatlie, C. T. Sachrajda, G. Martinelh, C. 
Pittori and G. C. Rossi, Nucl. Phys. B352, 
(1991) 266. 

2. P. Hasenfratz, Nucl. Phys. B (Proc. Suppl.) 
34, (1994) 3. 

3. R. C. Brower, unpubhshed CERN (1992). 

4. W. Bock, A. K. De, K. Jansen, J. Jersak, 
T. Neuhaus and J. Smit, Nucl. Phys. B344 
(1990) 207. 

5. R. C. Brower, Y. Shen and C-I Tan, "Chirally 
Extended Quantum Chromodynamics" , to be 
published in Phys. Rev. D (1994). 

6. S. Gottlieb, W. Liu, D. Toussaint, R.L. 



Renken and R.L. Sugar, Phys. Rev. D 35 
(1987) 2531. 

For a recent review, see, D. N. Petcher, Nucl. 
Phys. B (Proc. Suppl.) 30 (1993) 50. 



